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Abstract. Continuum observations at radio, millimetre, 
infrared and soft X-ray energies are presented for a sam- 
ple of 22 quasars, consisting of flat and steep spectrum 
radio loud, radio intermediate and radio quiet objects. 
The primary observational distinctions, among the differ- 
ent kinds of quasars in the radio and IR energy domains 
are studied using large observational datasets provided 
by ISOPHOT on board the Infrared Space Observatory, 
by the IRAM interferometer, by the sub-millimetre array 
SCUBA on JCMT, and by the European Southern Obser- 
vatory (ESO) facilities IRAC1 on the 2.2 m telescope and 
SEST. 

The spectral energy distributions of all quasars from 
radio to IR energies are analyzed and modeled with non- 
thermal and thermal spectral components. 

The dominant mechanism emitting in the far/mid-IR 
is thermal dust emission in all quasars, with the exception 
of flat spectrum radio loud quasars for which the presence 
of thermal IR emission remains rather uncertain, since 
it is difficult to separate it from the bright non-thermal 
component. The dust is predominantly heated by the op- 
tical/ultraviolet radiation emitted from the external com- 
ponents of the AGN. A starburst contributes to the IR 
emission at different levels, but always less than the AGN 
(< 27%). The distribution of temperatures, sizes, masses, 
and luminosities of the emitting dust are independent of 
the quasar type. 
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Radio continuum - In- 



1. Introduction 

Radio quiet and radio loud (not blazar) quasars (RQQ 
and RLQ, respectively) have similar spectral properties 
in the ultraviolet (UV), optical, and infrared (IR), but 
their radio powers differ by several orders of magnitude 



(Elvis ct al. 1994). This divergence takes place at millime- 
tre (mm) wavelengths. At these wavelengths the contri- 
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bution from two emission components merge, namely the 
synchrotron emission dominant in the radio domain and 
t hermal emission from cool d ust (30-50 K) in the far-IR 
( Barvainis fc Antonucci 1989|) . It is still not entirely clear 
whether the distinction between RLQ and RQQ is a conse- 
quence of differences in their central engines or whether it 
merely reflects differences in their environments. The pri- 
mary observational distinctions in the IR domain, and the 
proposed physical mechanisms to explain them are stud- 
ied here, using the new insights provided by Infrared Space 
Observatory^ (ISO; Kcsslcr ct al. 1996] ) measurements. 



1.1. The Radio 



emission 



Two main types of RLQ can be distinguished on the 
basis of their radio spectrum: the flat spectrum ra- 
dio loud quasars (FSRQ), and the steep spectrum ra- 
dio loud quasars (SSRQ). FSRQ show highly-collimated 
structures and very compact features, with flat or in- 
verted radio spectra. SSRQ have radio spectra dominated 
by synchrotron emission from extended radio lobes. The 
lobes and a radio core in the centre of these objects 
are signs of a relativistic jet. According to the unified 
scheme flBarthel 1989} |Urry fc Padovani 199$ FSRQ are 
the counterparts of SSRQ in which the jet is aimed at the 
observer. 

The origin of the much weaker radio emission in RQQ 
is far less certain. The majority of the total radio emis- 
sion from the RQQ comes from the compact features 
in the nucleus (< 1 kpc for unresolved regions, and at 
least 2 kpc for the resolved ones) rath er than the body 
of the host galaxy ( Kukula ct al. 1998| ). It has been pro- 
posed that the activity in RQQ is supplied by a star- 
burst, i.e. thermal bremsstrahlung and synchrotron emis- 
sion coming from strongly radiative supernova? and su- 
pernova? remnants (SNRs) in a very dense environment 



1 ISO is an ESA project with instruments funded by ESA 
Member States (especially the PI countries: France, Germany, 
the Netherlands and the United Kingdom) and with the par- 
ticipation of ISAS and NASA. 



2 



Polletta M. et al.: Far-IR of Radio Loud and Radio Quiet Quasars 



where shocks accelerate electrons ( Terlevich et al. 1992| ). 
Alternatively, if the energy supply arises from accre- 
tion onto a massive black hole, the radio emission from 
RQQ (as in RLQ) is caused by radio jets, but the 
bulk kinetic power of these jets is for some reason ~ 
10 3 times lower than those of RLQ flMiller et al. 199$ . 



This second hypothesis seems to be favored by recent 
studies, because of high brightness temperatures calcu- 
lated (typical SN/SNR s have T B < 10 5 K), by the evi- 
dence of a pc-scale jet ( |Blundcll fc Beasley 1998 ) and by 



observations of flat /in verted and variable radio spectra 
( Barvainis et al. 1996 ). 

Recently, quasars with intermediate radio luminosi- 
ties have been discovered and labeled Radio Intermediate 
Quasars (RIQ) ( [Francis et al. 1993| , [Falcke et al. 1995| ). 
RIQ may represent the Doppler boosted counter- 
parts of radio quiet quasars. This hypothesis is sug- 
gested by the variability observed at radio wavelengths 



(Falcke et al. 1996) 



1.2. The Infrared emission 

The presence of a dominant thermal (circumnuclear dust 
emission) , or non-thermal (synchrotron radiation from the 
AGN) component in the IR continuum of quasars is still 
debated. 

Many attempts to establish the origin of the 
IR emission in RLQ and RQQ have been per- 
formed through observations in the sub-millimetre 
(sub-mm) of quasars detected by IRAS (RLQ 



in |Chini et al. 1989a 


, and 


Antonucci et al. 1990; 


RQQ in Chini et al. 1989b, 


Barvainis et al. 1992, 


Hughes ct al. 1993, and Hughes et al. 1997; and both 


in Andrcani et al. 1999 


, this last work is the only one 



based on an optically selected sample). The main test 
applied to recognize the presence of thermal emission in 
the IR spectra of the objects was based on the slope of 
the continuum emission (F^ cx v a ) connecting the far-IR 
and sub-mm data. A steep, a > 2.5, continuum is strong- 
evidence for thermal dust emission. Most of the sources 
studied have a < 2.5 and are consistent with a dominant 
self-absorbed synchrotron emission component. However, 
some RQQ have spectral slopes as steep as a = 4.35, 
which, along w ith observation s of strong molecular gas 
(CO) emission ( Barvainis 1997 ), give strong support to a 
thermal mechanism as the origin of the far-IR component 
in RQQ. 

Among the RLQ, a is, at most, ~ 0.9 for the FSRQ, 
and 1.1 for the SSRQ flChini et al. 1989a|). Variability, 



shape of the continuum spectral energy distribution and, 
in some cases polarization, indicate that the radio, mm 
and far-IR emission of FSRQ is dominated by the syn- 



chrotron process (Lawrence et al. 1991). On the contrary, 
many SSRQ show evidence of thermal emission: their 
far-IR spectr a are brighter than e xtensions of the ra- 
dio emission ( Antonucci et al. 1990 ), suggesting a differ- 



ent origin than the non-thermal radio component; and the 
flux is constant, consistent with it arising from a region 



much larger than a light year (Edclson & Malkan 1987). 
Moreover, the spectral energy distributions of some RLQ, 
both FSRQ and SSRQ, show evidence for a galaxy 
component: reddening, residual starlight, molecular gas 
(Scoville et al. 1993), and some thermal dust emission 



in the near-IR ( Barvainis 1987| ). Both components, non- 
thermal synchrotron radiation and thermal dust emission, 
are probably present at IR wavelengths, as observ ed in the 
RLQ 3C273 (|Robson et al. 1986| ; |Barvaims 1987|) . 



1.3. Relation between the Radio and Infrared emission 

A tight, linear correlation is observed between 
t he far-IR flux and the radio fluxes in AGN 
( 5opp fc Alexander 1991 ), suggesting a common ori- 
gin. RQQ and RLQ occupy well defined regions in 
Log(i(IR))-Log(i(Radio)) space, and show a relation 
with a similar slope, just shifted to higher radio power 
by a factor ~10 4 in RLQ. RQQ show a similar relation 
as spirals, starbursts, and ultra luminous IR galaxies 
(ULIRG), suggesting that their IR emission may arise 
from sufficiently energetic star formation in the host 
galaxy ( {5opp fe Alexander 1991 ). However, the majority 
of the bolomctric luminosity in over half of known ULIRG 



seems to arise from a buried AGN (Sanders 1999). Addi- 
tionally, the variable and flat spectra, and high brightness 
t emperatures shown b y many RQQ at radio frequencies 
( Barvainis et al. 1996 ) suggest that the radio emission is 
related to the AGN rather than to a starburst. 



Proposed scenarios 



The unified model (Barthcl 1989; Urry fc Padovani 1995) 



predicts that similar disk-like dust distributions exist in 
both RQQ and RLQ. Orientation of the active nucleus, 
environment, and jet luminosity all affect the relative con- 
tributions of thermal and non-thermal sources to the ob- 



served infrared luminosity (Chini et al. 1989a) 



Other scenarios have been proposed to explain 
the large differences in radio power between RQQ 
and RLQ: different sp in of the central black hole 
( Wilson fc Colbert 1995 ), or different morphological type 
of the host galaxy. Indeed, different radio powers are ex- 
pected if one population of objects is fueled by mergers 
(ellipticals) and one is fueled by mostly internal processe s 
within the galaxy (spirals) ( Wilson fc Colbert 1995 ). 
However, recent studies on the host galaxies of quasars in- 
dicate that the host galaxies of RQQ are in several cases el- 
liptical and not always spiral galaxies (Taylor et al. 1996). 
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1.5. Open issues 

A better knowledge of the radio and IR properties of 
quasars is required to test the unified model predictions, 
and answer the following questions: 

1. What is the dominant mechanism emitting at IR wave- 
lengths in RLQ and RQQ? 

2. Do RLQ and RQQ have the same dust properties (tem- 
perature, source size, mass, and luminosity)? 

3. Does an interplay between the radio and the IR com- 
ponents exist? 

These questions can be addressed through the study of 
the spectral energy distributions (SED) of RLQ and RQQ. 
Here, we present the SEDs from radio to IR frequencies 
of a sample of 22 AGN (7 RQQ, 11 RLQ, 2 radio galaxies 
(RG) and 2 RIQ). The selected sample, even if incom- 
plete and heterogeneous, is useful to address these ques- 
tions thanks to several properties characterizing the sam- 
ple (steep/flat radio spectra, radio loudness/quietness), 
and to the large amount of photometric data available 
in the radio, mm/sub-mm and IR domains. This work is 
based mainly on IR data provided by ISO. ISO data re- 
duce the frequency gap between sub-mm and far-IR obser- 
vations, better sample the IR spectral band with a larger 
number of filters than previous instruments, and increase 
the number of detected objects thanks to a higher sensi- 
tivity. The study of the IR emission of quasars will be ex- 
tended in the future with the results of the European and 
of the U.S. ISO Key Quasar Programs providing a simi- 
lar coverage of the IR SED for a larger sample of quasars 
(see first results in Haas et al. (1998), and Wilkes et al. 
(1999)). 



Table 1. Basic properties of the sample 



Source RA (J2000) Dec (J2000) Type 1 z 



Name 


h 


m 


s 














3C 47 1 ' 4 


01 


36 


24 


+20 


57 


27 


SSRQ 





425 


PKS 0135-247 1 


01 


37 


38 


-24 


30 


53 


FSRQ 


o 


831 


PKS 0408-65 1 ' 2 


04 


08 


20 


-65 


45 


09 


RG 






PKS 0637-75 1 ' 2 


06 


35 


47 


—75 


16 


17 


FSRQ 


o 


653 


PG 1004+ 130 2 


10 


07 


26 


+ 12 


48 


56 


SSRQ 


o 


240 


PG 1048-090 2 


10 


51 


30 


-09 


18 


10 


SSRQ 


o 


344 


4C 61. 20 1 


10 


52 


33 


+61 


35 


20 


SSRQ 


o 


422 


PG 1100+772 1,5 


11 


04 


14 


+76 


58 


58 


SSRQ 


o 


312 


PG 1103-006 2 


11 


06 


32 


-00 


52 


52 


SSRQ 


o 


423 


PG 1216+069 1,2 


12 


19 


21 


+06 


38 


38 


RIQ a 


o 


331 


PG 1352+183 1 


13 


54 


36 


+18 


05 


18 


RQQ 





152 


PG 1354+213 1 


13 


56 


33 


+21 


03 


51 


RQQ 





300 


PG 1435-067 1 ' 3 


14 


38 


16 


-06 


58 


21 


RQQ 





129 


PG 1519+226 1 


15 


21 


14 


+22 


27 


43 


RQQ 





137 


PG 1543+489 1 ' 5 


15 


45 


30 


+48 


46 


09 


RQQ 





400 


PG 1718+481 1 


17 


19 


38 


+48 


04 


13 


RIQ a 


1 


084 


B2 1721+34 1 


17 


23 


20 


+34 


17 


58 


SSRQ 





206 


HS 1946+7658 1 ' 4 


19 


44 


55 


+77 


05 


52 


RQQ 


3 


020 


3C 405 1 


19 


59 


28 


+40 


44 


01 


NLRG 





056 


B2 2201+31A 1 


22 


03 


15 


+31 


45 


38 


FSRQ 





295 


PG 2214+139 3 


22 


17 


12 


+ 14 


14 


21 


RQQ 





066 


PG 2308+098 1 ' 3 


23 


11 


18 


+ 10 


08 


15 


SSRQ 





433 



t SSRQ: Steep Spectrum Radio Loud Quasar; FSRQ: Flat 
Spectrum Radio Loud Quasar; RQQ: Radio Loud Quasar; 
RIQ: Radio Intermediate Quasar; NLRG: Narrow Line Radio 
Galaxy; RG: Radio Galaxy. 

1 Observed with ISOPHOT; 2 Observed with SEST; 3 Ob- 
served with IRAC1; 4 Observed with IRAM; 5 Observed with 
SCUBA. 

a Classification in Falcke et al. (1996). 



2. Observational dataset 

Source names, coordinates, and redshifts of the se- 
lected sample were taken from the NASA Extragalac- 
tic Database (NED)Q, and are listed in Table [j]. 
Infrared observations were obtained for 18 of the 
sources with the Imaging Photopolarimeter on ISO 
(ISOPHOT; |Lemke et al. 1996Q , and 3 were observed 
with IRAC1 on the 2.2 m ESO/MPE telescope. Mil- 
limetre and sub-mm data were obtained for 10 objects 
with the IRAM interferometer at Plateau dc Burc in 



France ( Guilloteau et al. 1992 ), the Sub-millimetre Com- 
mon User Bolometer Array (SCUBA; |Holland et al. 1999|) 
on the James Clerk Maxwell Telescope^] and the Swedish 

2 The NASA/IPAC Extragalactic Database (NED) is oper- 
ated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and 
Space Administration. 

3 The James Clerk Maxwell Telescope is operated by the 
Joint Astronomy Centre on behalf of the Particle Physics 
and Astronomy Research Council of the United Kingdom, the 
Netherlands Organization for Scientific Research, and the Na- 
tional Research Council of Canada. 



ESO Sub-mm Telescope (SEST) of the European South- 
ern Observatory (ESO) at La Silla. The instruments used 
for each object are indicated in a footnote to Table [l]. 

The observational details (observing date, wavelength, 
and measured flux density) of IRAC1, SEST, IRAM, and 
SCUBA observations are reported in Table 0. The data ob- 
tained with the SEST telescope were reduced to outside 
the atmosphere, corrected for the gain elevation charac- 
teristic of the telescope, and calibrated with Uranus. All 
IRAM observations were performed in compact configu- 
ration. All five antenna were used during most of the ob- 
servations, with the exceptions indicated in a footnote of 
Table f| The calibrators were 3C 454.3 for 3C 47, and 
22040+420 and 1928+738 during the two observations of 
HS 1946+7658. SCUBA observations yielded good results 
at 850 /im only due to marginal weather. The fluxes were 
calibrated using the canonical gain value of 220 Jy/V, 
since a calibration scan at 850/im was not done. The cho- 
sen value is relatively insensitive to the weather and should 
be good to within 20%. ISOPHOT observations and data 
reduction are described in the next section. These new 
observations were supplemented with literature data from 
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Table 2. IRAC1, SEST, IRAM and SCUBA results Table 3. List of references of selected published data 



Source Name Instrum. 



Obs. Date A 
(d m y) (/im) 



pt 
(mJy) 



Source Name 



References number 1 



QP AH 


TP A T\/T 


i i 

14 


U f 


yo 


i Qfin 
loUU 


1 7 n+o o 
1 / .UztZ.o 






14 


07 


98 1 ' 3 


3000 


SO 8+0 fi 


PKS 0408—65 


SEST 


30 


11 


95 


1300 


<30 


Pl^Q f\«1T TKO 


O-CjO 1 


oU 


1 1 
1 1 


■Jo 


i Qfin 
loUU 


TO A 1 ^7 












1300 


Q05+S1 


PC 1004+ ISO 


SEST 


30 


11 


95 


1300 


<39 


PG 1048-090 


SEST 


30 


11 


95 


1300 


<18 


PG 1100+772 


SCUBA 


17 


01 


99 


850 


6.8+2.4 


PG 1103-006 


SEST 


30 


11 


95 


1300 


<27 


PG 1216+069 


SEST 


30 


11 


95 


1300 


<45 


PG 1435-067 


IRAC1 


21 


06 


96 


3.7 


12.1+0.4 












3.7 


12.2+0.5 


PG 1543+489 


SCUBA 


17 


01 


99 


850 


<5.4 


HS 1946+7658 


IRAM 


20 


05 


98 4 


1300 


<10.4 






20 


05 


98 4 


3000 


<5.0 


PG 2214+139 


IRAC1 


21 


06 


96 


4.7 


44.2+4.6 


PG 2308+098 


IRAC1 


21 


06 


96 


3.7 


<17.56 



* Upper limits to the flux are given at the 3<r 
observed on July 17, and 22. 2 Only four antenna 
July 17. 3 Only three antenna were used on July 
on July 17, and 22. 4 Also observed on May 23. 



level. 1 Also 
were used on 
14, and four 



radio to near-IR. We have also collected data at soft X- 
ray energies. For reasons of homogeneity we collected only 
ROSAT data, available for most of the sources. From pub- 
lished soft X-ray spectra we derived the flux at 1 keV cor- 
rected for absorption. The observed absorption is always 
compatible with the galactic absorption. The list of refer- 
ences from which data were retrieved, for each object, is 
reported in Table ||. 

3. ISOPHOT observations and data reduction 

Photometric data at several (up to 11) wavelengths be- 
tween 3.6 and 200 /im were obtained for each object using 
the single-element PI and P2 detectors plus the two array 
cameras, C100 (3 pixels x 3 pixels) and C200 (2 pixels x 
2 pixels). Detector and observing parameters are listed in 
Table [I| Most of the observations (124 in total for 16 ob- 
jects) were performed in chopper mode, and the remaining 
(37 for 10 objects) by mapping the region surrounding the 
target (scans or rasters). 

In chopper mode the radiation beam is deflected from 
the source (on-source position) to adjacent fields on the 
sky (off-source position) several times in order to measure 
the background emission. Triangular (T) and rectangular 
(R) chopping modes were used. In the triangular chopper 
mode the background emission is measured in two differ- 
ent regions, while in the rectangular chopper mode it is 



3C 47 

PKS 0135-247 
PKS 0408-65 
PKS 0637-75 
PG 1004+130 
PG 1048-090 
4C 61.20 
PG 1100+772 
PG 1103-006 
PG 1216+069 
PG 1352+183 
PG 1354+213 
PG 1435-067 
PG 1519+226 
PG 1543+489 
PG 1718+481 
B2 1721+34 
HS 1946+7658 
3C 405 

B2 2201+31A 

PG 2214+139 
PG 2308+098 



1, 2, 3, 4, 5 
1, 3, 5, 6, 7 
1 

1, 2, 3, 5, 6, 8, 9 

1, 10, 11, 12, 13, 14, 15, 16 

1, 3, 11, 14, 15, 17 

1, 3, 18 

1, 2, 3, 12, 15, 17, 19, 20, 21 

1, 12, 14, 15, 22, 23 

2, 12, 14, 23, 24, 25, 26 
2, 14, 15, 17, 25 

2, 14, 15, 25 
2, 14, 15, 25 
2, 14, 15, 25 

12, 14, 15, 25, 27, 28, 29 
1, 2, 14, 23, 26, 30 
1, 13, 17, 31 
25, 32 
1, 33, 34 

1, 2, 7, 13, 16, 20, 31, 35, 36, 37 

38, 39, 40, 41 

1, 12, 15, 29, 30, 42 

1, 15, 30 



van Bemmel et al. 1998 



Gezari et al. 1997 



t 1: NED; 2 

4: 

6: fTornikoski et al. 1996 



8: [Veron-Cetty fc Veron 1998 



10: 
12: 
14: 
16: 
18: 
20: 
22: 
24: 
26: 
28: 
30: 
32: 
34: 
36: 
38: 
40: 
42: 



Wilkes et al. 1994 



Banders et al. 1989 



Neu gebauer et al. 1987 



Ennis et al . 1982; ; 

1991 



Reid et al. 



Chini et al. 1989a; 



siebert et al. 1998 



Blundell fc Beasley 1998 



Brinkmann et al. 1997 ; 



Kiihr et al. 1981 ; 



Steppe et al. 1993|: 



Tanner et al. 1996; 



11: Kapahi 1991 



Falcke et al. 1996 



Barvainis et al. 1996 



Wang et al. 1996| 



Kuhn et al. 199 



Haas et al. 1998 



Hoekstra et al. 1997 



Neugebauer et al. 198E 



Neugebauer ct al. 197S 



Hughes et al. 1995 



13: 
15: 
17: 
19: 
21: 
23: 
25: 
27: 
29: 
31: 
33: 
35: 
37: 
39: 
41: 



Lister et al. 199 



Miller et al. 1993 



Elvis et al. 1994 



Antonucci et al. 1990 



Lonsdale & Morison 1983 



Kellermann et al. 198!: ; 



Yuan et al. 1998 



1 



Andreani et al. 



Chini et al. 1989b 



1999| ; 



Schartel et al. 1996 



Robson et al. 1998 
7999f 



Bloom et al. 



Mitchell et al. 1994 



Ghosh et al. 1994 



Ibrasranta et al~1992| : 



Table 4. ISOPHOT detector properties and covered sky 
region in mapping observations 



Detector 


A 


Pixel 


Scan 


Raster 


Name 


( A™ ) 


size 


Coverage 






(") 


( " x " ) 


( " x " ) 


PI 


3.6, 4.8, 7.3, 12 




52x156 




P2 


25 








C100 


60, 80, 100 


43.5 


138x230 


230x230 


C200 


150, 170, 200 


89.4 


184x460 


276x460 
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measured in only one position. Observing dates, filters, 
apertures, exposure times, chopping mode, and measured 
fluxes are reported in Table | for each 1SOPHOT chopper 
observations. 

In mapping mode the telescope moves in a pattern 
around the source, providing more sky coverage than in 
the chopper mode (Tabic [I]). PI detector maps were per- 
formed with an aperture of 52" during all observations, 
except one (B2 220I+31A) during which the chosen aper- 
ture was 23". Observing dates, filters, exposure times, and 
measured fluxes are reported in Table p. More details on 



mapping mode are reported in section p. 2.1 



3.1. First steps of the data reduction: from ERD to AAP 
level 

The first part of the data reduction was performed us- 
ing version 8.1 of the PHT Interactive Analysis (PIA)Q 
tool (Gabriel 1997). We started the reduction with the raw 
data processed with version 8.7 of the Off-Line Processing 



(Laureijs et al. 1998). The raw data form a sequence of 
detector read-outs distributed in 2™ (n=2-6) sets of four 
response curves or ramps, as function of time (Edited Raw 
Data: ERD in Volts). 

Each set of four ramps represents a sky position. Each 
ramp is corrected for the non-linearity of the detector re- 
sponse, and for contamination of cosmic particle events 
(glitches). The removal of read-outs affected by glitches is 
carried out by applying two median filtering techniques: 
the single-threshold technique that uses a threshold of 4.5 
standard deviations (ct) for flagging bad read-outs and the 
two-threshold technique that uses a threshold of 3. Oct for 
flagging and I.Oct for re-accepting read-outs. After apply- 
ing the non-linearity correction and the deglitching to the 
ERD, a straight line is fitted to each ramp, in order to 
determine its slope or Signal per Ramp Data (SRD in 
Volt/s). 

In most of the cases the first 25 or 50% (1 or 2 ramps 
of 4) of the signals per chopper plateau at the SRD level 
are discarded to enable the detector response to stabilize 
at the level corresponding to the source flux density. The 
remaining data are further corrected for highly discrepant 
points (value at more than 3ct from the average signal) 
still contaminated by glitches, for the orbital dependent 
dark current, and for the signal dependence on the ramp 
integration time (reset time interval) to obtain an average 
Signal per Chopper Plateau (SCP in Volt/s). 

After applying flat-fielding correction using PIA val- 
ues, the SCP data are calibrated to obtain the Standard 
Processed Data (SPD in unit of Watts). Since the detec- 
tor response varies with time, it is determined at the time 
of the observation by measuring the flux emitted by two 
thermal Fine Calibration Sources (FCS1 and FCS2) on 



board. The FCS measurements are reduced in the same 
way as the scientific measurements up to this step. Data 
from FCS1 are used because they are the best calibrated. 
The FCS1 signal is checked in order to remove data with 
large uncertainties (this step is equivalent to computing 
the weighted mean of the FCS1 data). 

In the case of mapping observations, the FCS1 is 
observed twice, before and after the observation of the 
source. The photometric calibration we use is the average 
value of the two FCS1 measurements. 

After the flux calibration the AAP (Auto Analysis 
Product) data are obtained. They are a sequence of 2 n 
off- and on-source flux measurements (in Jy) each cor- 
responding to a sky position. The reduction from the 
AAP level to the final results is performed using our own 
IDL routines, and not following the standard pipeline. 
This procedure was also applied in the reduction of 
ISOPHOT chopper data of a sample of Seyfert galaxies 



(Polletta fc Courvoisicr 1999) 



4 PIA is a joint development by the ESA Astrophysics Divi- 
sion and the ISOPHOT consortium. 



3.2. From the AAP level to final results 

The last steps of the data reduction before determining the 
source flux are the background subtraction, the deletion 
of remaining highly discrepant points, and the correction 
for effects depending on chopper plateau time, vignetting 
(only for chopper observations) and point spread function. 

In the case of chopper observations with the C100 de- 
tector only the central pixel pointed on the source is con- 
sidered to derive the flux density, since the eight border 
pixels contain only a small fraction of the central point 
source and summing these values would hence only in- 
crease the noise. 



3.2.1. Background subtraction 

In chopper observations the background is measured at 
each off-source position. Since in some cases the instru- 
ments show long term drift effects, the background sig- 
nal is estimated near the time of each on-source measure- 
ment and subtracted. In the case of chopper observations 
the background estimates are obtained by computing the 
weighted mean of each pair of consecutive off-source mea- 
surements. The weights are computed from PIA statistical 
uncertainties. Since the sequence of chopper plateaux ends 
with an on-source position, we used the weighted mean of 
the two last on-source measurements, and the flux ob- 
served in the last off-source position to determine the last 
pair of on- and off-source values, for a total of 2" _1 flux 
values. 

Small maps of the regions immediately surrounding ten 
of the targets were constructed in one or both of the follow- 
ing ways (Fig. Q): multiple linear scans across the source 
and rastering the detector about the source. A scan with 
the PI and CI 00 detectors consisted of three steps of the 
telescope, with this sequence repeated three times. Only 
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Table 5. Details of ISOPHOT Observations performed in 
Chopper Mode 



A 


Apt. 


Exp. 


Chop. 




Uncert. 1 ^ 


(H 


(") 


00 


Mode 


(mjy) Stat. 


Syst. 


3C 47 (January 30, 1997) 


4.8 


23 


128 


R 


<64. 






12.8 


23 


128 


R 


69. 


10. 


21. 


20 


23 


128 


R 


<171. 






60 




256 


R 


170. 


52. 


51. 


100 




256 


R 


246. 


52. 


74. 




PKS 0408- 


-65 (June 


s, 05 1997) 






4.8 


23 


128 


R 


<26. 






12.8 


23 


128 


R 


<60. 






20 


23 


128 


R 


<84. 






60 




128 


R 


<308. 






100 




128 


R 


<162. 








PKS 0637- 


-75 (June 


s, 05 1997) 






4.8 


23 


128 


R 


<26. 






12.8 


23 


128 


R 


<43. 






20 


23 


128 


R 


<99. 






60 




128 


R 


<117. 






100 




128 


R 


<143. 






4C 61.20 (April, 27 1996) 


7.3 


13.8 


512 


T 


<35. 






12 


23 


256 


T 


<34. 






25 


52 


512 


T 


<110. 






60 




128 


R 


<278. 






80 




128 


R 


258. 


48. 


77. 


100 




128 


R 


<303. 






150 




128 


R 


<267. 






170 




128 


R 


<148. 






200 




128 


R 


<335. 







In units of mjy. 



the middle row of C100 pixels (8, 5, and 2, as depicted 
in Fig. |l|) viewed the source. The C200 scans contained 
four steps, repeated twice, with the source centered be- 
tween two pixels. Note that the source was observed in 
only the middle two steps of the C200 scan. The raster 
patterns were 3x3, 3x3, and 4x2 or 2x4, for the PI, 
C100, and C200 detectors, respectively. Each pixel viewed 
the source once in the raster maps. The step size between 
exposures for both scans and rasters was approximately 
equal to the pixel size, a little more to take into account 
the gap between pixels, for the C100 and C200 cameras, 
and equal to the aperture size for the PI detectors, which 
resulted in a different total sky coverage for each detector 
(see Table A background estimate for each on-source 
measurement was obtained from a weighted average of the 
flux measured in the raster or scan positions immediately 
preceding and following the source position by the same 
pixel. Using the same pixel to determine the background 



Table 5. (continued) 



A 


Apt 


. Exp. 


Chop. 


F„ 


Uncert.^ 


(Urn) 


(") 


00 


Mode 


(mjy) Stat. 


Syst. 




PG 1100+772 (June 


■., 17 1996) 






4.8 


52 


512 


R 


25. 


7. 


8. 


7.3 


52 


512 


R 


32. 


5. 


10. 


12 


52 


512 


R 


36. 


8. 


11. 


25 


120 


512 


R 


54. 


10. 


16. 


60 




256 


R 


70. 


14. 


21. 


100 




128 


R 


<282. 






150 




512 


R 


<417. 






200 




512 


R 


<468. 






PG 1216+069 (July, 11 1996) 


4.8 


52 


512 


R 


50. 


15. 


15. 


7.3 


52 


512 


R 


72. 


6. 


22. 


12 


52 


512 


R 


90. 


12. 


27. 


25 


120 


512 


R 


95. 


26. 


29. 


60 




128 


R 


65. 


14. 


20. 


100 




128 


R 


<140. 






150 




512 


R 


<137. 






200 




512 


R 


<149. 






PG 1352+183 (December, 14 1996) 


60 




64 


R 


198. 


23. 


59. 


80 




64 


R 


<237. 






100 




64 


R 


<213. 






150 




64 


R 


<283. 






170 




64 


R 


<240. 






200 




64 


R 


<206. 








PG 1354+213 (June 


, 13 1996) 






7.3 


13.8 


i 2048 


T 


18. 


5. 


5. 


12 


23 


1024 


T 


21. 


5. 


6. 


25 


52 


2048 


T 


23. 


7. 


7. 


60 




128 


R 


<298. 






80 




128 


R 


<194. 






100 




128 


R 


<125. 






150 




128 


R 


<170. 






170 




128 


R 


<137. 






200 




128 


R 


<280. 







In units of mjy. 



reduced the impact of uncertainties in the flat field. The 
weighted average background was subtracted from each 
on-source measurement providing a sequence of source flux 
values. 

Residual effects of detector instabilities, in both chop- 
per observations and maps, produced occasional dis- 
crepant points, which were culled by one-pass sigma clip- 
ping. The threshold number of standard deviations to re- 
ject a flux value depended on the number of points in the 
sequence (Chauvenet's criterion in Taylor (1982)), ranging 
from 1.15cr to 2.66er. 
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Table 5. (continued) 



A 


Apt 


Exp. 


Chop. 




Uncert.* 


(um) 


(") 


(8) 


Mode 


(mJy) 


Stat. 


Syst. 




PG 1435- 


067 (January, 07 1997) 




7.3 


13.8 


256 


T 


70. 


9. 


21. 


12 


23 


64 


T 


<94. 






25 


52 


128 


T 


<718. 






60 




64 


R 


<1111. 






80 




64 


R 


<849. 






100 




64 


R 


<333. 






150 




64 


R 


<386. 






170 




64 


R 


<369. 






200 




64 


R 


<317. 






PG 1519+226 (February, 01 1997) 


7.3 


13.8 


256 


T 


32. 


8. 


10. 


12 


23 


64 


T 


<48. 






25 


52 


64 


T 


76. 


17. 


23. 


60 




128 


R 


172. 


44. 


52. 


80 




128 


R 


<194. 






100 




128 


R 


121. 


22. 


36. 


150 




128 


R 


<207. 






170 




128 


R 


<151. 






200 




128 


R 


<154. 






PG 1543+489 (May, 30 1996) 


4.8 


52 


512 


R 


20. 


6. 


6. 


7.3 


52 


256 


R 


39. 


4. 


12. 


12 


52 


256 


R 


45. 


6. 


14. 


25 


120 


512 


R 


140. 


13. 


42. 


60 




128 


R 


470. 


144. 


141. 


100 




128 


R 


399. 


55. 


120. 


150 




128 


R 


455. 


117. 


182. 


200 




256 


R 


<377. 






PG 1718+481 (May, 30 1996) 


4.8 


52 


512 


R 


18. 


3. 


5. 


7.3 


52 


256 


R 


24. 


6. 


7. 


12 


52 


256 


R 


25. 


6. 


8. 


25 


120 


512 


R 


44. 


7. 


13. 


60 




128 


R 


90. 


1. 


27. 


100 




128 


R 


56. 


1. 


17. 


150 




128 


R 


<251. 






200 




256 


R 


<261. 







In units of mJy. 
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Table 5. (continued) 



A 


Apt. 


Exp. 


Chop. 


F„ 


Uncert.* 


(Mm) 


(") 


(s) 


Mode 


(mJy) Stat. 


Syst. 


B2 1721+34 (April, 20 1996) 


7.3 


13.8 


512 


T 


19.* 


5. 


6. 


12 


23 


512 


T 


<26.* 






25 


52 


512 


T 


<135.* 






HS 1946+7658 (May, 27 1996) 


4.8 


52 


1024 


R 


<5. 






7.3 


52 


1024 


R 


4.5 


1.2 


1.4 


12 


52 


1024 


R 


8.7 


0.7 


2.6 


25 


120 


1024 


R 


47. 


15. 


14. 


60 




64 


R 


<230. 






100 




64 


R 


<300. 






150 




64 


R 


682.* 


17. 


273. 


200 




64 


R 


378.* 


27. 


151. 


B2 2201+31A (November, 16 1996) 


3.6 


5 


2048 


T 


<17. 






4.8 


7.6 


1024 


T 


<52. 






7.3 


13.8 


256 


T 


45. 


5. 


14. 


12 


23 


64 


T 


64. 


8. 


19. 


12 


23 


64 


T 


<62. 






25 


52 


128 


T 


106. 


25. 


32. 


60 




128 


R 


193. 


13. 


58. 


80 




128 


R 


139. 


11. 


42. 


100 




128 


R 


212. 


51. 


64. 


150 




128 


R 


847. 


97. 


339. 


170 




128 


R 


699. 


98. 


280. 


200 




128 


R 


382. 


97. 


153. 


PG 2308+098 (November, 26 1996) 


3.6 


5 


2048 


T 


<31. 






4.8 


7.6 


1024 


T 


13. 


2. 


4. 


12 


23 


128 


T 


<39. 






PG 2308+098 (December, 12 1996) 


7.3 


13.8 


256 


T 


<21. 






12 


23 


128 


T 


16. 


55. 


5. 


25 


52 


512 


T 


<108. 






60 




128 


R 


<183. 






80 




128 


R 


<174. 






100 




128 


R 


<83. 






150 




128 


R 


<201. 






170 




128 


R 


<152. 






200 




128 


R 


<214. 







In units of mJy. * Doubtful data. Discussed in Sec.|3.3|and 4.1 



In the case of chopper observations with the C200 cam- 
era, the source flux was computed by adding together the 
fluxes measured by each of the four pixels. The source 
flux was divided between pairs of adjoining pixels in the 
C200 scans; these were averaged by weighting with their 
uncertainties. Table ^ lists the weighted mean of each 
flux sequence after clipping, or 3cr upper limits for non- 
detection, where a is the quadratic sum of the statisti- 



cal and systematic uncertainties (see section pT3| ) , also re- 
ported in the tables. 

3.2.2. Vignetting correction 

In the case of chopper observations with the C100 and 
C200 detectors, the data are further corrected for the sig- 
nal loss outside the beam of the telescope (vignetting). 
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Table 6. Details of ISOPHOT Observations performed in 
Raster Mode 



Obs. date 


A Exp. 


F„ 


Uncert. 1 " 




H ( s ) 


(mJy) 


Stat. 


Syst. 


3C 47 


January, 31 1998 


60 636 


<129. 






PKS 0135-247 


December, 12 1997 


12 d 1206 


<199. 






PKS 0408-65 


June, 05 1997 


170" 195 


<51. 






PKS 0637-75 


June, 05 1997 


170" 346 


<180. 






March, 15 1998 


150" 387 


<171. 








200° 387 


<203. 






PG 1100+772 


August, 15 1997 


60 348 


<99. 








100 348 


<55. 








150 b,c 579 


<82. 








200 b ' c 579 


<54. 






October, 28 1997 


12 d 954 


<161. 








60 d 630 


68. 


22. 


20. 


November, 02 1997 


150" 579 


<120. 






PG 1543+489 


November, 01 1997 


12 d 864 


<202. 








60 d 630 


269. 


27. 


81. 




150 a ' d 574 


279. 


45. 


84. 



* In units of mJy. "4x2 map. ft 2x4 map. c Step amplitude of 
90". d Telescope nodding mode. 



The PIA default vignetting correction factors were com- 
puted considering the dependency only on the distance of 
the chopper positions and on the filter, but recent inves- 
tigations have shown that they depend also on the time 
per chopper plateau (M. Haas, private communication). 

The PIA default vignetting factors were applied di- 
rectly to the data, resulting in little change to the flux 
values obtained with the C100 camera. However, the fluxes 
from the C200 detector are ~ 80%, in median, larger after 
the vignetting correction. Given the uncertainty over the 
accuracy of the vignetting factors, we averaged the fluxes 
derived with and without the correction and added the 
respective uncertainties in quadrature. 

3.2.3. Correction for effects depending on the chopper 
plateau time and point spread function 

In the case of chopper observations a correction for the sig- 
nal loss due to effects depending on the time per chopper 
plateau is applied to the computed flux value. Since short 
integration times with the C100 and C200 detectors do 
not reach the full signal, the observed flux can be reduced 



Table 6. (continued) 



Obs. date 


A 


Exp. 


F„ 


Uncert.^ 




w 


(s) 


(mJy) 


Stat. 


Syst. 


PG 1718+481 


April, 20 1997 


60 


348 


59. 


16. 


18. 




100 


348 


40. 


13. 


12. 




150" 


323 


<104. 








200 a 


579 


<56. 






November, 01 1997 


12 d 


864 


<164. 








60 d 


630 


69. 


13. 


21. 




150 a,d 


574 


<145. 






HS 1946+7658 


April, 28 1997 


60 


348 


49. 


12. 


15. 




100 


348 


79. 


14. 


24. 




150" 


451 


<114. 








200 a 


579 


<356. 






October, 26 1997 


200° 


579 


<227. 






November, 17 1997 


12 d 


954 


<171. 






November, 18 1997 


12 d 


954 


<168. 








60 d 


630 


66. 


15. 


20. 


December, 08 1997 


60 d 


630 


68. 


20. 


20. 


3C 405 


October, 30 1997 


60 d 


630 


2360. 


76. 


708. 




l7Qa ,d 


574 


1582. 


160. 


475. 


B2 2201+31A 


November, 20 1997 


12 d 


852 


55. 


15. 


17. 




60 d 


414 


111. 


20. 


33. 




170" 


387 


<249. 







f In units of mJy. "4x2 map. b 2x4 map. c Step amplitude of 
90". d Telescope nodding mode. 



by large factors, typically up to 68 % for the C100 detec- 
tor, and up to 12% for the C200 detector for the shortest 
observations. The correction factors we use are the PIA 
default values. 

All the computed fluxes are finally corrected for the 
point spread function (psf) of each detector using the de- 
fault PIA values derived empirically in most of the cases. 
The available psf correction values correspond to a source 
position centered in a pixel or located in a corner of the 
pixel. In case of scans with the C200 detector we needed 
the psf correction corresponding to the target located in 
the middle of a side of the pixel. We derived it by assuming 
a bi-dimensional Gaussian function for the psf and con- 
straining its parameters using the other two known values 
for each wavelength (30.6% at 150 /an, and 29.1% at 170 
fj,m). 

3.3. Calculation of systematic uncertainties 

During the data reduction only statistical errors were 
taken into account, and not those in the absolute flux den- 
sity calibration. The accuracy of the absolute photometric 
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Telescope Nodding Classical Raster 

P1-P2 Detector (3 x 3) 



C100 Detector (3 x 3) 



C200 Detector (4 x 2) 



Fig. 1. Schematic representation of the two observing 
modes used for mapping observations: telescope nodding 
(scan) on the left side and classical raster on the right 
side. The position of the source in the different cases is 
represented by a star. The bold square represents the ini- 
tial position of the detector. The numbering corresponds 
to the pixels numbers. The detectors move first from the 
left to the right, then in the opposite direction. In the 
case of classical raster, at each change of the horizontal 
direction, the detector shifts down one step. All the steps 
are separated by the pixel size. The reported pixel sizes 
do not correspond to the real proportions. The numbers 
in parenthesis refer to the number of horizontal steps x 
the number of repetitions for nodding and vertical steps 
of telescope motion in the raster maps. 

calibration depends mainly on systematic errors (detector 
transient effects, calibration response, dark current, point 
spread function), and it is currently known to be better 
than 30% (Klaas et al. 1998). The associated statistical 
and systematic uncertainties are reported in the last two 
columns of Table [| and ^. We associate to the measured 
flux an uncertainty that is the quadratic sum of their sta- 
tistical uncertainties and 30% (for the C200 measurements 
obtained in chopper mode we use 40% that corresponds to 
half of the uncertainty due to the vignetting correction) of 
the measured value. Among our ISOPHOT observations 
the datasets of B2 1721+34 were identified as failed after 
inspection of their quality by the ISO team because they 



were heavily affected by cosmic rays and thus scientifically 
useless. We report these data in Table ||, and indicate that 
they are doubtful in a footnote. 

4. The spectral energy distributions 

We report the first far-IR detection for 9 (4C 61.20, PG 
1216+069, PG 1352+183, PG 1354+213, PG 1435-067, 
PG 1519+226, PG 1718+481, PG 2308+098, and HS 
1946+7658) of the 18 sources observed with ISOPHOT. 
Among the remaining 9 sources, three (PKS 0135—247, 
PKS 0408-65, and PKS 0637-75) were not detected, and 
the remainings 6 were already detected by IRAS. 

All the data were converted to monochromatic lumi- 
nosities in the rest frame of the object (Ho = 75 km s _1 
Mpc -1 , qo=0.5) using the following equation: 



= ^diF vohs /{\ 



(1) 



where -F^obs is the monochromatic flux in the observer's 
frame, and c?l is the luminosity distance to the object. 
In the case of PKS 0408—65 we adopted a redshift value 
equal to 0.5, arbitrarily chosen, since no redshift measure- 
ment is available. This source will not be considered in the 
following analysis. 

The spectral energy distributions (SEDs) as vh v ver- 
sus v in the rest frame of all sources are shown in Fig. ^[ 
Upper limits are plotted only if there is no detection 
at that frequency at any epoch, and in cases of mul- 
tiple upper limits at the same frequency from different 
epochs, only the most stringent is used. A broad band 
spectrum of a local galaxy in its rest frame, represent- 
ing the AGN host galaxy, is superposed in each panel 
of Fig. ||. We chose the spiral galaxy M100, modeled by 
Silva et al. (1998), to which we added data in the ra- 
dio ( Becker et al. 1991 ; |Gregory~ 



Condon 1991; White 



Becker 1992) and in the soft X-ray energy domains 



(Immlcr ct al. 1998) in the case of RQQ, and the template 
of a giant elliptical galaxy (Silva et al. 1998) in the case 
of RLQ, RIQ and RG. The reported host galaxy template 
was not modified by any normalization. In many cases 
this is orders of magnitude below the observed luminosi- 
ties, and even if it is shifted towards higher luminosities 
to reach the quasar SED, it will remain below at most 
frequencies. 



4--1- The nature of the IR emission: thermal or 
non-thermal emission ? 

The IR emission can have a thermal or non-thermal origin. 
Several investigation methods can be applied to identify 
the emission process: 

1. The value of the slope of the sub-mm/far-IR spectral 
break can discriminate between optically-thick (self- 
absorbed) synchrotron and thermal emission from dust 
grains. The self-absorbed synchrotron model is char- 
acterized by a maximum value of the sub-mm/far-IR 
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10 L 



PKS0135-247 



PKS0637-752 (High) 



B22201+31A 



1989 



PKS0637-752 



PKS0637-752 (Low) 




— I 1 1 — l h 

P22201+31A (Low) 




10" 



10 



10 



10 11 10 
v (Hz) 



12 



10 



13 



10 14 10 E 



10" 



10 



10 



10 11 10 
v (Hz) 



12 



10 



13 



10 



14 



Fig. 2. SEDs as vL v versus v in the rest frame of the objects. All sources shown on this page are FSRQ. Full circles 
represent literature data, stars SEST data, open diamonds ISOPHOT chopper data, triangles and filled diamonds 
ISOPHOT raster data at different epochs, and squares IRAS data. Arrows represent 3<r upper limits. The solid line 
represents the sum of all single fitted components repre sent ed by dotted lines. The dotted lines represent the best fit 
non-therma l mo dels of the radio component (see section L2, and Table 0), and the parabolic fits of the IR component 
(see section 5.2). When more than one radio component is present a symbol is reported for each spectral component: 
C for core and L for lobe (Li and L2 for two lobes). A typical host galaxy template (dashed line), a spiral galaxy in 
the case of RQQ, and a giant elliptical galaxy in the case of RLQ, RIQ and RG, is over-plotted. The radio spectrum 
of PKS 0637—752 was fit in two ways, using only the simultaneous mm data corresponding to the flattest power law 
fit (c) (dash-dot line in b) and the two corresponding to the steepest power law fit (d) (dash-dot-dot-dot line in b). 
The radio spectrum of B2 220f +3iA was also fit in two ways, taking all data from radio to mid-lR (e) and only data 
from radio to far-IR (f). The flattest power law fit (dash-dot line in e) and the steepest one (dash-dot-dot-dot in e) 
measured during simultaneous mm observations of B2 2201+31A are reported. The big open circles in e correspond 
to simultaneous observations at mm and near-lR wavelengths. 



slope that is 2.5 (F„ tx v a ), if the radiation is emit- 
ted by an electron population with a simple power- 
law energy distribution, or somewhat larger (a < 3) 
if a thermal electron pool or dual power-law energy 
distribution is invoked ( ic Kool fc Bcgclman 1989| : 
[Bchlickciscr ct al. f99f ). In cither case the maximum 
synchrotron slope is attained only for a completely 
homogeneous source, otherwise it will be lower. The 



asymptotic thermal sub-mm/far-IR slope is expected 
to be > 3 since the optically thin thermal spectrum 
derives from Rayleigh- Jeans law with an additional 
parameter dependent on frequency, F v cc B U (T) ■ 
where r<j, the dust optical depth, is oc with (3 ~l-2 



(see section 4.3). 

A non-thermal origin of the IR emission is indicated if 
relatively short time scale flux variability is observed, 
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Fig. 2. Continued. SEDs for SSRQ are displayed on this page. 



since a synchrotron component is expected to come 
from a very compact source. Most of the dust emission 
comes from an extended source with a long variability 
time scale. 

A thermal origin can be attributed to the IR emission 
if a non- varying excess from any reasonable extrapola- 
tion from the radio domain in the plot vL v versus v is 
observed (Hughes et al. 1997). 



It is possible to distinguish the origin of the IR emission 
based on brightness temperature (Sanders et al. 1989) 
and polarization measurements, but the lack of these kind 
of data for the sample does not allow us to apply them. 

The first test will prove the thermal origin of the far- 
IR emission only if the coldest dust component is also the 
brightest one. Colder and less bright dust components will 
flatten the sub-mm/far-IR spectral slope. Two sources in 
the sample, PG 1543+489 and 3C 405, have a sub-mm/far- 
IR spectral index larger than 2.5. Their far-IR emission is 
therefore dominated by thermal radiation. The remain- 
ing sources have insufficient data at long wavelengths to 
constrain the sub-mm/far-IR spectral slope. 



The observation of no variability is not conclusive, 
while if a flux variation is observed the non-thermal hy- 
pothesis will be strongly supported. This method can be 
applied only to those sources that were observed several 
times. In our sample, at least two IR observations (IRAS 
and ISOPHOT, or several ISOPHOT observations) at dif- 
ferent epochs are available for ten sources (3C 47, PKS 
0637-752, PG 1100+772, PG 1543+489, PG 1718+481, 
B2 1721+481, HS 1946+7658, 3C 405, B2 2201+31A, and 
PG 2214+139) (see Fig. g). The data relative to the same 
observation epoch, instrument and observation mode are 
represented with the same symbol and connected by a line 
in Fig. g. At least two measurements at the same wave- 
length are available for all the ten objects with the ex- 
ception of PG 2214+139. Six sources (PG 1718+481, PG 
1100+772, 3C 405, PG 1543+489, B2 1721+34, and PG 
2214+139) show no sign of variability. Two of them also 
satisfy the first test. For two other sources (3C 47, and 
B2 2201+31A) the observed variation is only marginally 
<1.6er, where a includes the statistical and the system- 
atic uncertainties. We consider their emission as constant 
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0) 



in the span of our observations. In the case of PKS0637- 
752 we can only give a lower limit of the variation since 
the source was not detected by ISOPHOT. At wavelengths 
shorter than 60^m ISOPHOT and IRAS give consistent 
results, while at longer wavelengths they differ of more 
than 1.9cr. 

The 150/im flux values for the high redshift source HS 
1946+7658 differ by at least a factor of 6 between two 
epochs separated by less than a year (see Tables || and ^) . 
Such extreme FIR variability is unlikely in a RQQ. The 
larger flux resulted from a chopped measurement, which 
is more susceptible at longer wavelengths to both back- 
ground fluctuations and instrumental effects (e.g., a num- 
ber of inadequately corrected cosmic ray strikes, unusually 
high detector drift, etc..) than the raster maps. Fig. [| 
shows a l.lxl.l degree region around HS 1946+7658 at 
100/im from IRAS, with the position of the source and 
the ISOPHOT chopper direction indicated. Bright cirrus 
structures are nearby, and the background ranges from 5 IRSKY is An Observation Planning Tool for the Infrared 
0.21 - 0.38 MJy/sr within 90"of the source (equivalent to Sky developed for NASA at IPAC, JPL/Caltech. 
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Fig. 2. Continued. 



0.16 - 0.28 Jy on the C200 covered area ) (IRSKY0 ver- 
sion 2.5). A strong gradient in the direction SE-NW in 
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the background emission is clearly shown, which coincides 
with the ISOPHOT chopper direction. These variations 
may be greater at longer wavelengths and with the better 
spatial resolution of ISOPHOT. As a consequence of the 
uncertainties for this source, chopped data at A > 100/im 
for HS 1946+7658 will not be included in subsequent anal- 
ysis. 

The result of this method suggests that for all the se- 
lected objects, with the exception of PKS 0637—752, the 
IR emission is thermal in origin. The result is not very 
strong since a non-thermal process may also produce a 
constant flux. However, it is unlikely to measure the same 
flux from a non-thermal source during three different ob- 
servations performed in a range of 14 years (from 1983 to 
1997) as observed in six objects (see Fig. ||). 

The third test is the simplest, and it can be applied 
to the whole sample since a broad spectral coverage from 
radio to the IR is available for all the sources. Before ap- 
plying this test we need to estimate the contribution from 
the non-thermal component in the IR and subtract it from 



the observed IR spectrum. The non-thermal contribution 
can be estimated by fitting the radio data with a reason- 
able model and then extrapolating it to higher frequencies 
(see next section). 

4-2. Contribution of the radio non-thermal component in 
the infrared 

In order to estimate the contribution of the radio non- 
thermal component in the IR domain, the radio continuum 
was fitted with some plausible models, and extrapolated 
to higher frequencies. In the following we will distinguish 
two components in the radio spectra of RLQ, the extended 
component (radio lobes), and the core component (unre- 
solved core, jet, etc.). 

The radiation emitted by compact sources, like the 
cores and the hot-spots of radio loud objects, can be mod- 
eled by a self-absorbed synchrotron emission spectrum. In 
the case of a homogeneous plasma with isotropic pitch an- 
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Fig. 4. IRAS observation at 100 /im of a sky region of 
l.lxl.l degrees around the source HS 1946+7658, whose 
position is indicated by S. The sky region where the back- 
ground was measured is indicated by B. The white poly- 
gon represents approximatively the region observed by the 
C200 camera. The solid line indicates the chopper direc- 
tion -41° to the West of North. 

gle distribution and power law energy distribution of the 
form N(E) oci? _s , this can be expressed as: 



L v oc [ — 



1 - exp ( - ( - 



(2) 



where u t is the frequency at which the optical depth of 
the plasma is equal to unity, and f cu toff is the frequency 
corresponding to the cut off energy of the plasma energy 
distribution, at which the energy gains and losses of the 
electrons are equal. The optically thick and optically thin 
spectral indices are denoted by ai, and ct2- In the case of 
a homogeneous source, a\ is expected to be 2.5, and «2 is 
related to the exponent s of the plasma energy distribution 
by the relationship: «2 = — (s — l)/2. The superposition 
of several self-absorbed components can produce a flatter 
power law. In this case the spectral model of equation (|J) 
will remain valid, but ai will not have the same meaning. 
If the source is optically thin, as in extended sources (ra- 
dio lobes), the emitted spectrum can be expressed more 
simply as: 



L„ oc 



fcut off 



. e "cut off 



(3) 



In many cases the contribution of the synchrotron com- 
ponent at high frequencies is negligible compared to the 
observed emission, therefore a high energy cut off was not 
included in the model. The model is then represented by 
a broken power law of slopes a\ and a%, or by a simple 
power law of slope 0:2, according to the presence/ absence 
of self absorption. 



The synchrotron model describes well the observed ra- 
dio spectra of most of the sources in the sample. However, 
while in the case of SSRQ it is quite easy to separate the 
different components and hence apply a model for each of 
them, the spectral modeling is more difficult for FSRQ. 
For these sources we parameterize the emitted spectrum 
with an empirical equation that is valid if the resulting 
spectrum is produced by self-absorbed synchrotron emis- 
sion or by optically thin synchrotron emission due to a 
hard electron spectrum produced through the acceleration 
processes in turbulent plasma ( Wang et al. 1997 ). The ob- 
served fiat radio spectra are described by the following 
equation: 



L v oc 



1 — exp 






1 — exp 







(4) 



where v\, is the frequency at which the spectrum flat- 
tens, a is the spectral index observed at low frequencies 
(y < z/t,)i an( i the other model parameters have the same 
meaning as in equation (^). 

The observed SEDs from the radio to the mm energy 
domains were fitted with one or a combination of these 
models (see best fit parameters in Table |?]) . In some cases 
near-IR data have also been used in the fits, in particular 
when no data in between mm and near-IR frequencies con- 
strained the spectrum to lie below the near-IR flux (PKS 
0135-247, PKS 0637-752, 4C 61.20, PG 1004+130, PG 
1216+069, PG 1718+481, 3C 405, and B2 2201+31A). 
In a few cases we fixed some model parameters, as in- 
dicated in a footnote of table [?], since the available data 
could not constrain them. The fixed values were chosen 
in the range of values that provided reasonable spectra 
with properties similar to those observed in other sources. 
Model a in Table corresponds to equation (||). It was 
applied for modeling core spectra of PG 1216+069, and 
3C 405. The same model without the cut off (6) was ap- 
plied for modeling weak core spectra for which the high 
energy cut off was not constrained (3C 47, PG 1004+130, 
4C 61.20, PG 1048-090, PG 1100+772, PG 1103-006, 
and PG 1718+481). A simple power law model (c) was 
used to fit the radio emission from the lobes of SSRQ and 
RG (3C 47, PKS 0408-65, PG 1004+130, 4C 61.20, PG 
1048-090, PG 1100+772, PG 1103-006, B2 1721+34, 3C 
405 and PG 2308+098) and the radio emission of RQQ 
(PG 1543+489, and PG 2214+139). Simultaneous obser- 
vations available in the literature generally do not provide 
wide or well-sampled wavelength coverage, so all avail- 
able data were used in the fits to the SEDs. The data 
and analysis are adequate for the central purpose of es- 
timating the contribution of the non-thermal component 
to the IR emission. The model d, corresponding to equa- 
tion (Q), was applied to fit the radio emission of FSRQ 
(PKS 0135-247, PKS 0637-752, and B2 2201+31A). The 
value of the break frequency Vh was arbitrarily fixed to 
2.75 GHz, since it provides a good fit to the emitted spec- 
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Table 7. Best fit parameters of non-thermal models 



Source Name 
Source 


ivioaei 


ota 


Oil 


Oil 


Vt V cut off 

(10 9 Hz) (10 14 Hz) 


3C 47 (Core) 


b 




1.64 


-0 


45 


2.98 




3C 47 (Lobe) 


c 






-0 


99 






PKS 0135-247 


d 


-0.78 


0.30 


-0 


80 


69.52 


2.00 


PKS 0408-65 


c 






-1 


19 






PKS 0637-752 


d 


-0.60 


0.55 


-0 


78 


37.20 


4.49 


PKS 0637-752 (Low) 


e 


-0.55 


0.55 


— 1 


50 


75.12 




PKS 0637-752 (High) 


d 


-0.60 


0.49 


-0 


78 


41.28 


4.22 


PG 1004+130 (Core) 


b 




0.53 


-0 


19 


2.00 




PC 1004+130 (Lobe) 


c 






-0 


80 






4C 61.20 (Core) 


b 




0.20 


-0 


15 


20.00 (F) 




4C 61.20 (Lobe 1) 


c 






-1 


04 






4C 61.20 (Lobe 2) 


c 






-0 


93 






PG 1048-090 (Core) 


b 




0.50 


-0 


65 


20.00 (F) 




PG 1048-090 (Lobe) 


c 






-0 


84 






PG 1100+772 (Core) 


b 




0.15 


-0 


80 


11.00 




PG 1100+772 (Lobe) 


c 






-0 


83 






PG 1103-006 (Core) 


b 




0.30 


-0 


•50 


20.00 (F) 




PG 1103-006 (Lobe) 


c 






-0 


70 






PG 1216+069 


a 




0.76 


-0 


09 


0.96 


20.39 


PG 1543+489 


c 






-0 


90 






PG 1718+481 


b 




0.58 


-0 


34 


11.45 




B2 1721+34 


c 






-0 


69 






3C 405 (Core) 


a 




0.50 (F) 


-0 


74 


126.00 (F) 


3.58 


3C 405 (Lobe 1) 


c 






-1 


26 






3C 405 (Lobe 2) 


c 






-1 


30 






B2 2201+31A (Low) 


d 


-0.52 


0.60 


-0 


30 


28.60 


0.034 


B2 2201+31A 


(I 


-0.44 


0.50 


-0 


50 


34.89 


0.24 


PG 2214+139 


c 






-0 


92 






PG 2308+098 


c 






-0 


85 







* Model a corresponds to equation (|2j) , and model b to the same equation without cut off; 
model c corresponds to a simple power law; model d to equation (W), and model e to the 
same equation without cut off. 
(F) indicates a fixed value. 



trum of the three objects. In the case of PKS 0637—752 
(Low) (see section 4.2.1) the cut off was not included in 
the model (e) since the high frequency part of the spec- 
trum is very steep. 



4.2.1. Uncertainties in the radio contribution estimate 

The location of the high energy cut off is difficult to estab- 
lish. Every power law relative to the optically thin emis- 
sion was extended at higher frequencies until the spec- 
trum turned down, and hence a cut off was required by 
the data. A spectral cut off was thus required only in five 
objects (PKS 0135-247, PKS 0637-75, PG 1216+069, 3C 
405, and B2 2201+31A), but it could have been located at 
lower frequencies and present in other objects, too. In most 
of the cases this parameter does not affect the presence 
and the strength of the remaining IR flux, but its energy 
value may be important in FSRQ (PKS 0135-247, PKS 
0637-75, and B2 2201+31 A), since these objects have flat 



radio spectra for which extrapolation up to IR frequencies 
is comparable to the IR fluxes. For these sources a more 
accurate analysis of their radio spectra is needed. Since 
PKS 0135—247 was not detected in the IR, no further 
analysis can be performed. We concentrate only on PKS 
0637-75 and B2 2201+31A. In order to better constrain 
the non-thermal radio spectrum, i.e. to find some evidence 
of a spectral cut off at sub-mm/far-IR frequencies, we 
searched in the literature for simultaneous observations 
at these wavelengths, and we selected those that showed 
the flattest and the steepest spectrum. For PKS 0637—75 
the flattest mm po wer law, chosen among several simulta- 
neous observations ( Tornikoski ct al. 1996 ) , was measured 
on February 15th, 1990 (a(3.0-1.3 mm) = -0.77), and the 
steepest one was measured on April 4th, 1991 (a(3. 0-1.3 
mm) = —1.47). The two power laws are reported in Fig. ||b 
with a dashed, and a dashed-dotted line, respectively, plus 
displayed separately with flattest (Fig. ||c) and steepest 
(Fig ||d) spectral fits. The flattest spectrum overlaps the 
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observed IR spectrum, leaving no additional IR compo- 
nent. On the contrary, the extrapolation of the steepest 
spectrum to IR frequencies is clearly below the observed 
IR spectrum, but the IR observations were not simultane- 
ous to the mm observations. The source was observed by 
IRAS in 1983, and by ISO at different wavelengths in 1997. 
During the elapsed time the source became fainter in the 
far-IR, while shorter wavelength data from the two difer- 
ent epochs are consistent. In the following we will suppose 
that a thermal IR component is present, but dominating 
only at A <60 /im, and we will analyze its properties and 
compare them with those observed in other sources. 

For B2 2201+31A the flattest mm power law (a(1.0- 
0.87 mm) = —0.0 9) was measured on February 1989 
( |Chini ct al. 1989"a ), and the steepest one was measured 
on September 14th, 1993 (a(2.0-l. 3-1.1 mm) = -0.72). 
The two power laws are reported in Fig. ||e with a dashed, 
and a dashed-dotted line, respectively. The spectrum is 
in both cases quite flat, however the extrapolation of the 
1993 spectrum lies below the IR spectrum. More than the 
sub- mm data, the analysis of the emission at shorter wave- 
lengths gives important indications on the origin of the IR 
emission. B2 2201+31A was observed on September 15th, 
1993 also in the near-IR (simultaneous sub-mm and near- 
IR data are indicated by large open circles in Fig. ^Js). 
The near-IR data are above the extrapolation of the sub- 
mm data, suggesting the presence of two different spec- 
tral components in these two wavelength ranges (see the 
analogous case of 3C 273 in Robson et al. (1986)). This 
hypothesis is also suggested by the constant emission ob- 
served up to 60 /mi. A non-thermal source is expected 
to vary more at higher frequencies, due to greater energy 
losses. All these considerations suggest the short wave- 
length continuum is dominated by a thermal component. 
As in the case of PKS 0637—75 we will suppose that an 
additional IR thermal component is present at A < 60 
/im. 

These two sources (PKS 0637-752, and B2 2201+31A) 
are good examples of how variability can create an artifi- 
cial IR spectral turnover, or hide a real one. An IR spec- 
tral turnover may be due to different luminosity states of 
the source at different epochs, instead of to the presence 
of a separate IR component. The weakness of the radio 
emission in RQQ precludes that its extrapolation could 
account for the IR emission for all reasonable assump- 
tions on the radio variability. In SSRQ the extrapolation 
of the radio component in the IR is usually too faint to 
explain the IR emission, even if we take into account vari- 
ability. The variability factors observed in two SSRQ in 
our sample, 3C 47 and PG 1004+130, are too small to ex- 
plain the much higher IR fluxes, and this is probably true 
for the SSRQ in general. In the mm domain we measured 
a flux variation from the core of the SSRQ 3C 47 of a 
factor of ~ 2 in almost three years (the emitted flux den- 
sity at ~ 100 GHz was equal to 16 .3+0.9 mJy on Septem- 
ber 1995 ( |van Bemmel et al. 1991) , and equal to 30.8±0.6 



mJy on July 1998 (this work)). The SSRQ PG 1004+130 
was observed twice at 6 cm, in 1982 and in 1984, with a 
flux variation of a factor ~2.5, from 12 mJy to 30 mJy 



(Lister et al. 1994) 



In conclusion, the radio models shown in Fig. ^| indi- 
cate the presence of an additional IR component in almost 
the whole sample. According to the third test, this result 
indicates that the observed IR emission is of thermal ori- 
gin. The properties of the IR e missi on in quasars will be 
derived and analyzed in section L3, after subtraction of 
the non-thermal contribution extrapolated from the radio 
domain. 



4-3. Modeling of the IR component 

The IR emission can be accounted for by reradiation of the 
central luminosity by gas and dust in warped disc s in the 



host galaxies of the quasars ( Sanders et al. 1989 ), in the 



outer edge of the accretion disc and in a torus of molecular 



(Niemeyer & Biermann 1993 




Granato & Danese 1994; 


Granato et al. 1997 




Pier & Krolik 1992], 1993), and/or 


by starburst emission 


( Rowan- Robinson 1995]). The host 



galaxy starlight contribution is probably negligible in the 
far/mid-IR since the host galaxy spectrum largely differs 
in shape and luminosity from the SED of the selected 
objects (see Fig. ||). We describe here the main obser- 
vational properties of the different objects of each class 
and compare them using a very simple model of thermal 
emission: the grey body model. This model does not 
take into account the source geometry (toroidal, warped 
disc, etc). An isothermal grey body at the temperature 
T emits at frequency v a luminosity density given by the 
following equation ( pear 19'88| , Wccdman 1986): 



47rV-B(2/ em) T)-(l 



(5) 



where r is the radius of the projected source, B(z/,T) is the 
Planck function for a blackbody of temperature T, and Td 
is the optical depth of the dust. The optical depth can be 
approximated by a power law of type t<j = (f/fo) t where 
Uq is the frequency at which the dust becomes optically 
thin, and j3 is the dust emissivity index. A non-linear least 
squares fit was used in the fitting procedure, leaving the 
radius r, the temperature T, and the frequency Uq free to 
vary , while the emissivity exponent [3 was fixed equal to 



1.87 ( [Polletta fc Courvoisier 1999| ). 

The observed IR SEDs are smooth and indicate a wide 
and probably continuous range of dust temperatures, de- 
scribable by several grey body components. The best fit 
grey body models of the observed IR SEDs are shown 
in Fig. |^. The thick solid line represents the sum of non- 
thermal and grey body components. Each individual com- 
ponent is represented by a dotted line. The temperature 
(T) and the size (r) of each grey body component are 
listed in columns 5-10 of Table @. It is worth noting that 
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Fig. 5. SEDs as vL v versus v in the rest frame of the objects (2 FSRQ, 2 RIQ, 1 RG, 1 RQQ). Symbols as in Fig. |[ 
Dotted lines represent the best fit non-thermal models of the radio component and the best fit grey body models of 
the IR component. The temperature of each grey body component is reported. The sum of all single fitted models is 
shown by a solid line. 



Table 8. Best fit parameters of grey body models 



Source Name Total L(IR) L ^L(z-^pan) Log(Md) I Component II Component III Component 
(10 11 L a ) (Ma) T (K) r (pc) T (K) r (pc) T (K) r (pc) 



3C 47 


36.1 


0.205 


7.55 


42.7 


7621 


208 


109 






PKS 0637-752 (Low) 


58.6 


0.056 


6.68 


61.0 


1923 


219 


276 






PC 1004+130 


11.6 


0.040 


5.35 


92.6 


663 


421 


12 






4C 61.20 


29.8 


0.086 


6.14 


86.8 


1401 










PC 1100+772 


14.1 


0.024 


5.38 


74.9 


597 


174 


49 


696 


2.16 


PC 1103-006 


25.3 


0.024 


5.18 


120.0 


496 






1322 


0.33 


PC 1216+069 


28.8 


0.007 


4.76 


115.8 


310 


511 


8 






PC 1352+183 


3.8 


0.091 


5.26 


80.9 


556 






1584 


0.11 


PC 1354+213 


9.3 


0.001 








211 


23 


898 


0.98 


PC 1435-067 


2.8 












4 


465 




PC 1519+226 


3.8 


0.097 


5.97 


50.7 


1264 


168 


41 


591 


1.34 


PC 1543+489 


49.9 


0.202 


7.54 


46.0 


8997 


165 


183 


561 


5.17 


PC 1718+481 


105.3 


0.007 


5.06 


141.5 


573 






878 


3.60 


B2 1721+34 


8.6 


0.047 


5.36 


80.0 


611 


230 


32 


1800 


0.11 


HS 1946+7658 


761.5 


0.006 


6.00 


130.8 


1294 


528 


37 






3C 405 


4.6 


0.137 


5.79 


67.7 


1265 


184 


13 






B2 2201+31A (Low) 


19.5 


0.002 








222 


18 


611 


3.99 


PC 2214+139 


2.0 


0.091 


5.22 


63.7 


540 


399 
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1900 


0.06 
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Fig. 5. Continued (6 SSRQ). 
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we could fit the observed IR spectra using a different opti- 
cal depth function (different /3, and Uq values) . The optical 
depth value is important in a discussion of the source ge- 
ometry in terms of an extended or compact heating source. 
In our models the optical depth values derived by the fits 
are low (< 1) in the far/mid-IR, and ~1 in the near-IR (~ 
3 /im). If the dust becomes optically thin at longer wave- 
lengths, the real source sizes will be smaller than our esti- 
mates, and vice versa. Using our optical depth values, the 
estimated sizes of the observed dust components range be- 
tween 0.06 pc and 9.0 kpc, and the temperatures between 
43 K and 1900 K. The minimum temperature may be due 
to an absence of dust at large distances (few kpcs) or at 
low temperature, and/or to starlight heating to the or- 
ders of the inferred minima. The maximal temperature 
is generally explained as a drop in opacity caused by the 
sublimation of the most refractory grains at temperatures 
T - 2000 K ( [Sanders et al. 1989|) . The total luminosi- 
ties observed in the IR, obtained by integrating the grey 
body components (see column 1 in Table ||), vary over a 
wide range, from 2.0- 10 11 L Q to 7.6- 10 13 L G . No signifi- 
cant difference in the distribution of sizes, temperatures, 
and luminosities are observed among different types of 
quasars. We also derive the mass of each dust component 



at the measured temperature, using the following equation 



(Hughes et al. 1997) 



M d = 



K d { 

k'emis) -B(z- / em i s , T)' 

where ifd(^emis) oc v@, /3 = 1.87, is the rest-frequency dust 
absorption coefficient. The normalization is Kd(v em i s ) = 
10 cm 2 g~ 1 at 250 /mi ( |Hildebrand 1983[ ), giving K A (v) = 
1.14 cm 2 g _1 at 800 /im. The range of assumed values 
of K<\(v) a t 800 /im in the literature is 0.4-3.0 cm 2 g _1 
([Draine fc Lee 1984 [Mathis fc Whiffen 1989j). Our dust 



mass estimates can thus differ by, at most, a factor 2.7. 
The derived values of dust masses are reported in column 
4 of Table ||, and, separately for each class, in Fig. ^. Since 
the largest dust masses are located in the outer, less illumi- 
nated, lower temperature regions of the dust distribution, 
Md is mainly constrained by far-IR data. Therefore, when 
sub-mm and far-IR data are not available, the real dust 
mass cannot be well measured. For this reason we did not 
report the dust and gas masses when the low tempera- 
ture component was not constrained. The absence of data 
in the near-IR has a negligible effect on the dust mass 
estimate. As for the other parameters (T, L(IR), r), the 
dust mass distribution does not differ significantly among 
different types of quasars (see Fig. H). 
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Fig. 5. Continued (6 RQQ). 



,13 



10 J 
v (Hz) 



10 



14 



10 



15 




Log( M, ) ( M. ) 

Fig. 6. Histogram of the dust masses for the different 
classes: RQQ, SSRQ, FSRQ, and RIQ. 



5. Similarities and differences in the SED of RLQ 
and RQQ 

The present sample contains a range of radio source clas- 
sifications, with which we can elucidate the dependence of 
broad-band spectral features on radio properties, thereby 



testing some unification scenario predictions. The limita- 
tions of these tests lie in the sample's relatively small size 
and heterogeneous nature. 

5.1. Average SED 

A quick look at the main spectral differences between the 
different kinds of quasars is provided by the comparison 
of the average SED of RQQ, RIQ, FSRQ, and SSRQ (in- 
cluding the RG 3C 405). The SEDs are shown in L u -u and 
v\j v -v spaces separately for each class over the radio/soft 
X-ray frequency range in Fig. |^. The broad spectra of two 
typical host galaxies (a giant elliptical and a spiral galaxy), 
in their rest frames and without any normalization, are 
also plotted in Fig. ^, as in Fig. ||. The average SEDs have 
been computed using the conventional mean, excluding 
upper limits. The width of each frequency bin is equal to 
0.5 in Log(z/). The reported uncertainties correspond to 
the standard deviation of the mean of the data per fre- 
quency bin. All the data have been connected by straight 
lines. At soft X-ray energies we indicate the average power 
law ± ltr computed from the distribution of best fit soft 
X-ray power law models of all objects of the same class 
(photon index T = 2.73±0.61 (RQQ), 2.39±0.19 (RIQ), 
2.39±0.23 (SSRQ), and 2.25±0.12 (FSRQ)). 
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Fig. 7. The average quasar energy 
distribution for RQQ (solid line), 
RIQ (dotted line), FSRQ (dashed 
line), and SSRQ (dash-dot line). 
The number of objects that were 
used in each class is reported on the 
figure with the class name. A sym- 
bol is reported for each frequency 
bin, if at least one data point is 
available. The two high redshift ob- 
jects of the sample: HS 1946+7658, 
and PG 1718+481, were excluded 
since they are characterized by very 
high luminosities that modify the 
average SEDs at a few frequen- 
cies producing an irregular spectral 
shape. Two host galaxy templates 
are also reported: a typical spiral 
galaxy (dash-dot-dot-dot line), and 
a giant elliptical galaxy (long dashed 
line). No normalization was applied 
to the reported curves. 



As expected, the largest difference in luminosity among 
the different classes appears at radio wavelengths. A 
smaller difference is observed at soft X-ray energies, and 
in the near-IR (y > 10 14 Hz corresponding to A < 3 /una), 
while the luminosity and the spectral shape in the mid- 
and far-IR are remarkably similar (see Fig.0). The large 
difference in the IR spectral shape between quasars and 
the host galaxy templates indicates that the contribution 
from the host galaxy is negligible also at radio and soft 
X-ray energies, and not only in the far/mid-IR. This re- 
sult is in agree ment with previous s tudies on the br oad 
SED of quasars QSandcrs et al. 1989j ; |Elvis et al. 1994| ). A 
quantitative comparison of the luminosity emitted at dif- 
ferent frequencies by each quasar class is presented in the 
next section. 



5.2. Multi band luminosities 

The IR component was also modeled by fitting a parabola 
in Log L^-Log v space (see Fig. ||). The parabola model 
gives a rough estimate of the strength and shape of the IR 
component, even if the spectral coverage is not complete. 
For several objects upper limits were also used in the fit. 
This model has by itself no physical meaning, however, it 
describes the IR component relatively well, it can easily 
be traced even with poor spectral coverage, and can take 
into account the whole IR emission of most of the objects 
in a larger wavelength range than the detailed grey body 
models. The parabola is too narrow to satisfy the observed 
IR SED in a few cases, e.g., in 3C 405 and PG 1543+489. 
In these cases we fitted only the far / mid-IR data where the 
IR emission usually peaks. The parabola parameters are 



its width, the frequency of maximum luminosity density 
(^peak) and the maximum luminosity density (L Upesik ). The 
parabola fit to the IR component was applied to all objects 
of the sample, except PKS 0135—247, for which no IR data 
are available (see Fig. |^a). The distribution of the peak 
frequencies values (fpeak) observed in the four different 
classes of objects is reported in Fig. || The distribution is 
quite similar for SSRQ and RQQ, ranging from 2.6T0 12 Hz 
(114 (Mm) to 3.6- 10 13 Hz (8 fim), while it is shifted to higher 
frequencies for FSRQ and RIQ, ranging from 9.0T0 12 Hz 
(33 (Mm) to 2.8-10 13 Hz (11 /im). This difference may be 
due to the flat radio non-thermal component extending to 
high frequencies in FSRQ and RIQ, and dominating the 
dust emission. 

We define the IR luminosity L(IR) as the product of 
the luminosity value at which each parabola peaks and 
the corresponding frequency (L(IR) = ^ pC ak • ^i/ pcak ) ■ Note 
that this parameter does not depend on the width of the 
parabola. Only upper limits for L(IR) could be derived 
for PKS 0408-65 and PG 1040-090. The distribution of 
L(IR) is reported in Fig. |[ In this, and in the following 
histograms upper limits are shown with arrows, one per 
object. The similarity in the IR luminosities and spectra 
(see also Fig. 0) in all quasars suggest a similar origin. 

The radio emission in the RLQ arises from two very 
different spatial scales, the core and extended components. 
We calculated the average of vh v over the rest-frame in- 
terval 5-9 GHz for each spatial component in all of the 
RLQ, except PG 1354+213 and HS 1946+7658, which 



were undetected at these frequencies. Fig. 10 displays his- 
tograms for the two components separately, and Fig. [ll] 
shows the distribution of the median of all measured vh,. 
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Fig. 8. Histogram of the IR peak frequency of the 
parabola model for the different classes: RQQ, SSRQ (in- 
cluding the RG 3C 405), FSRQ, and RIQ. 




Log( L(IR) ) ( erg s"' ) 



Fig. 9. Histogram of the peak luminosity of the IR 
parabola model for the different classes: RQQ , SSRQ 
(including the RG 3C 405 whose position is indicated), 
FSRQ, and RIQ. 



over the same frequency range, without component dis- 
tinction. The distribution of the median radio luminosity 
is bi- modal (Fig. |ll|). However, if we consider only core 
radio luminosities (top panel of Fig. |l^), the SSRQ radio 
luminosity distribution shifts towards lower values, mak- 
ing a continuous distribution, rather than a bi-modal one, 
but without overlapping. The contribution from the ex- 
tended components are very similar in FSRQ and SSRQ 
(bottom panel of Fig. |l0|) . In the following analysis we will 
consider only the core luminosity L (Radio). When the core 
luminosity is not available (PKS 0408-65, B2 1721+34, 
3C 405, and PG 2308+098), we report an upper limit cor- 
responding to the average radio luminosity relative to the 
extended component. 
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Fig. 10. Histogram of the core and extended components 
radio luminosity for the different classes: RQQ, SSRQ 
(including the RG 3C 405 whose position is indicated), 
FSRQ, and RIQ. 
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Fig. 11. Histogram of the median radio luminosity for the 
different classes: RQQ, SSRQ (including the RG 3C 405 
whose position is indicated), FSRQ, and RIQ. 

In the soft X-ray, we define L(SX) as vL v with v cor- 
responding to 1 keV in the observer's rest-frame. The dis- 
tribution of L(SX) for each class is reported in Fig. [l2| 
In the soft X-ray, no data are available for 3C 405, and 
PKS 0408—65, and only an upper limit is available for PG 
1004+130. 

5.3. Origin of the observed luminosities 

The main factors determining the observed luminosities 
are: the energy emitted by the central engine (AGN); the 
amplification due to Doppler boosting in a relativistic jet; 
and the contribution from a starburst. We will estimate 
the role of each of these parameters in producing the SEDs 
through the comparison of the observed radio, IR, and soft 
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Fig. 12. Histogram of the soft X-ray luminosity at 1 keV 
for the different classes: RQQ, SSRQ, FSRQ, and RIQ. 

X-ray luminosities, represented by L (Radio), L(IR) and 
L(SX), respectively (see section 5.2 for their definition). 



5.3.1. Orientation effects in RLQ 

The orientation of the beamed emission can be estimated 
from the radio core fraction R. This quantity, defined as 
the ratio between the core radio luminosity and the lumi- 
nosity of the extended radio emission at 5-9 GHz in the 
rest frame, serves as an orientation indicator of the radio 
source with respect to the obser ver, measuring t he relative 
strength of the core component ( Hes et al. 1995 ). The core 
flux was not available for three SSRQ and one RG, thus 
the parameter R was not computed. In the case of FSRQ 
we computed the luminosity of the extended component in 
the frequency range 5-9 GHz by extrapolating the power 
law observed at low frequencies (power law index ao in 
Table 0). 

The FSRQ are well separated from SSRQ in the dis- 
tribution of the ratio R (see Fig. |l3|). This difference per- 
mits us to estimate the enhancement factor of the beamed 
emission after a few considerations. First, the observed 
radio emission in RLQ is mainly produced by the jet 
and its core rather than a starburst, since star-emitting 
ULIRG have much lower radio luminosity than RLQ 
flColina fc Pcrcz-Olea 199§ ). Second, we assume that the 
radio source is intrinsically identical in FSRQ and SSRQ, 
and that the difference in their radio emission is due only 
to the orientation of the beamed emission. After these ap- 
proximations we can write 



L(Radio Core) FSR Q = A ■ L(Radio Core) SSR Q 



(7) 



where A is the amplification factor of the beamed emis- 
sion. Since the luminosity of the extended components are 
the same for the fiat and steep radio quasars (see above), 
using equation (Q) we can derive a relation between the 
parameters -Rfsrq and -Rssrq: -Rfsrq = A ■ i?ssRQ- Re- 
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Fig. 13. Luminosity emitted in the soft X-ray, radio and 
IR domains compared to the parameter R given by the 
ratio between the radio luminosity of the core and of the 
extended components. The different symbols correspond 
to different classes of objects: squares for SSRQ, and cir- 
cles for FSRQ. The dashed lines are the best fit lines. 

Table 9. Correlation test results for FSRQ and SSRQ 



X 


Y 


N 


n 


Pi (%) 


Log(L(IR)) 


Lo; 


5(R) 


7 


0.53 


25. 


Log(L(Radio)) 


Lo; 


;(R) 


9 


0.94 


0.03 


Log(L(SX)) 


Lo; 


;(R) 


8 


0.73 


6.0 



placing the observed values of i?ssRQ (— 0.05-0.15), and 
-Rfsrq (— 3-4) (see Fig. |l3|) in the above relation yields 
A -20-80. 

Fig. [ll] displays i(Radio), L(SX), and L(IR) versus 
the core fraction R. Linear correlation results for these 
relations are reported in Table [)[ where the parameter 
pairs are reported in the first two columns, the number of 
data pairs in the third, in column 4 the linear correlation 
rank (rj), and in column 5 the associated probability to 
have such a correlation rank from uncorrelated values (Pi). 

Higher radio and soft X-ray luminosities are observed 
in objects with higher values of the radio core fraction 
R (when the jet points towards us). The orientation ef- 
fect is more important in the radio domain, as shown by 
the stronger correlation, than in the soft X-ray, and neg- 
ligible in the IR. This implies that the radio core and a 
fraction of the total emitted soft X-ray luminosities are 
emitted anisotropically. We furthermore verified that R is 
not correlated with the redshift and thus that the above 
result is not an artifact of distance related biases in the 
measurement of R. 
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Assuming that the soft X-ray source is intrinsically 
identical in FSRQ and SSRQ, the observed difference in 
L(SX) arises from the orientation of the fraction, /, that is 
beamed. If the fraction of emitted radiation that is beamed 
is enhanced by a factor A, identical to that of the radio 
emission, the following relation between the soft X-ray 
luminosity in FSRQ and that in SSRQ will be valid: 

^(SX) FSRQ = A- f ■ L(SX) ssm + (1 - /) ■ L(SX) 

SSRQ' (8) 

Using the average value of the ratio 
L(SX) FgRQ /L(SX) SSRQ that is w 3.2 ± 2.2, and the 
range of values obtained for the factor A, we derive a 
fraction / « 3-12% for the beamed fraction of the soft 
X-ray component. 



FSRQ 




Log( L(60-200/xm)/L(3-6(Vm) ) 



5.3.2. SSRQ versus RQQ 

The radio and the soft X-ray luminosities are mainly pro- 
duced by the AGN component (see sections [5j] and 5.3.1 ) . 
The comparison between the luminosities emitted in the 
radio and soft X-ray domains is then equivalent to a com- 
parison of the AGN power in the two types of quasars. The 
radio core emission of SSRQ is on average 200 times higher 
than that of RQQ (see Fig. 10) and the soft X-ray lumi- 
nosity is on average 8 times higher in SSRQ than in RQQ 
(see Fig. |l2|). Since the SSRQ show luminosities higher 
than RQQ not only at radio and soft X-ray energies, but 



also in the hard X-ray domain (Lawson & Turner 1997), 
we argue that the bolometric AGN luminosity is much 
higher in SSRQ than in RQQ. The difference in the AGN 
power should be observable at all frequencies where the 
AGN emission dominates. We have already pointed out 
the similarity in IR luminosities and spectra of SSRQ and 
RQQ (see Figs. and ||). This similarity suggests that the 
origin of the dominant IR component is not AGN-related. 
The candidate is then a star burst. 

Some indication of the dominant IR emission mecha- 
nism can be gleaned from the shape of the SED. An AGN 
can emit a significant fraction, often a majority, of its in- 
frared luminosity at shorter wavelengths, A < 60/im, as 
long as the obscuring columns are not so large as to be 
optically thick at these wavelengths. Starburst dominated 
galaxies, on the other hand, produce the bulk of their in- 
frared emission at A > 60/im. The ratio of the luminosi- 
ties in these two wavelength regimes, L(60-200/im)/L(3- 
60/im), thus provides a rough estimate of the primary 
driver of the infrared component. A histogram of this 
ratio is presented in Fig. |l4], in which only the sources 
having at least two grey body components with T < 
1000 K are included. For comparison, the luminosity ratio 
from an average SED of low reddening starburst galax- 



ies (Schmitt et al. 1997) is also indicated. All of the AGN 



in the sample have infrared luminosity ratios less than the 
starburst fiducial value (=0.76) by a factor or four or more 
(the maximum ratio is 0.20 corresponding to 27% of star- 
burst contribution), suggesting that the infrared in these 



Fig. 14. Histogram of the ratio L(60-200/mi)/L(3-60/mi) 
for the different classes: RQQ, SSRQ, FSRQ, and RIQ. 
The dashed line represents the value derived from the av- 
erage SED of a sample of low reddening starburst galaxies. 

sources is dominated by the central engine. The RQQ and 
SSRQ have similar average ratios, 8% of the total IR emis- 
sion is produced by a starburst. 

Since the dominant IR source is the AGN, which dif- 
fers in SSRQ and RQQ, other factors have to be taken into 
account to explain the similarity of their IR luminosities 
and average spectra. These factors can be the dust cover- 
ing factor, its amount and geometric distribution. A higher 
dust covering factor in RQQ, due to larger total dust mass 
and/or a particular geometric distribution, may account 
for the similar IR radiation. Aspects of this hypothesis are 
attractive, since the presence of more circumnuclear mate- 
rial may be related to the physical conditions that hinder 
the jet formation and/or development in RQQ. However, if 
larger dust covering factors were present in RQQ, a higher 
probability to observe a RQQ through optically thick dust 
would be expected and a large fraction of RQQ with ab- 
sorption features in the soft X-ray would be observed. Ob- 
servations do not support these predictions (see section ||). 
Moreover, the similarity in the observed IR luminosities 
and spectra requires a fine tuning between the dust heat- 
ing source, the dust amount and its geometric distribution, 
which makes this explanation improbable. We propose an- 
other scenario in which the dust properties (amount and 
distribution) and the heating source are similar in SSRQ 
and RQQ. The dust distribution contributing in the far- 
IR probably extends from the more external regions of 
the AGN, and is predominantly heated by the optical and 
UV radiation fields filling these external regions which are 
similar in both classes. At relatively large distances from 
the centre, the AGN components are similar in both types 
of quasars, towards the innermost regions, where the soft 
X-rays are emitted and a jet is formed, SSRQ and RQQ 
become different. The high energy photons escape from 
the centre without significant dust absorption and provide 
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an important probe of the central radiation source. Ob- 
servations indicate that the soft X-ray radiation is higher 
in SSRQ (see Fig. ^) in agreement with the proposed sce- 
nario. 

From this analysis it is suggested that the main differ- 
ence between RQQ, and SSRQ takes place in the inner- 
most nuclear regions where the emitted power is higher in 
SSRQ than in RQQ, while the AGN external regions (dust 
distribution and optical/UV source) have similar proper- 
ties in both types of quasars. 

Since FSRQ and SSRQ show similar properties in the 
IR, once the non-thermal contribution is subtracted, this 
conclusion can be also extended to FSRQ. 

6. Conclusions 

Continuum observations from radio to soft X-ray energies 
have been presented for a sample of 22 quasars character- 
ized by different degrees of radio emission. The IR data 
were obtained with ISOPHOT and IRAC1, and the mm 
data with IRAM, SEST and SCUBA. Further IR and mm 
data, and some radio and soft X-ray fluxes, were drawn 
from the literature. The availability of a broad band SED 
for several types of quasars allows us to separate the dom- 
inant spectral components in the radio and IR energy 
bands, and thus to compare their spectral properties in 
different types of quasars. The spectral analysis and the 
comparison of luminosities emitted in the radio, IR, and 
soft X-ray energy bands yields the following results: 

1. What is the dominant mechanisms emitting at IR en- 
ergies in RLQ and RQQ? 

In our quasar sample the dominant mechanism emit- 
ting in the far/mid-IR is thermal emission from dust 
heated by the optical and UV radiation produced in 
the outer regions of the AGN. A starburst contributes 
to the IR emission at different levels, but always less 
than the AGN (<27%). The presence of thermal IR 
emission in FSRQ remains rather uncertain. Among 
the three FSRQ of the sample we cannot derive any 
conclusion for two of them, and for a third one (B2 
2201+31A) the observed data suggest a dominant ther- 
mal component at A < 60 fim. 

2. Do RLQ and RQQ have the same dust properties (tem- 
perature, source size, mass, and luminosity)? 

The equilibrium temperature of dust grains, the size 
and the mass of the dust distribution, and the emitted 
luminosity have been evaluated for all quasars. The 
estimated sizes of the observed dust components lie 
between 0.06 pc and 9.0 kpc, and the temperatures be- 
tween 43 K and 1900 K. The total luminosity observed 
in the IR, obtained by integrating the grey body com- 
ponents, varies over a wide range: 2-760T0 11 L Q . The 
amount of emitting dust in all types of quasars also 
varies in a broad range : 6-10 4 -4-10 7 M Q . The distribu- 
tion of any of these parameters does not differ signifi- 
cantly among the different types of quasars. 



3. Does an interplay between the radio and the IR com- 
ponents exist? 

A bright and flat non-thermal component can be suf- 
ficiently strong in the IR to mask the dust emission in 
some sources, particularly FSRQ. However, this does 
not mean that the dust emission is absent. After sub- 
tracting an IR synchrotron component extrapolated 
from the radio, the residual IR emission had similar 
spectral shape and luminosity, regardless of the radio 
properties. 

These results are based on the analysis of a small sam- 
ple, and have to be confirmed by the study of larger sam- 
ples. ISOPHOT has doubled the number of quasars with 
IR detections in the sample presented here. A great deal of 
additional progress on understanding the IR properties of 
quasars is expected when all of the quasar data available 
in the ISO archive will be fully analyzed and studied. 
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